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Abstract. We study the dynamics of ELKO in the context of accelerated phase of our 
universe. To avoid the fine tuning problem associated with the initial conditions, it is required 
that the dynamical equations lead to an early-time attractor. In the earlier works, it was 
shown that the dynamical equations containing ELKO fields do not lead to early-time stable 
fixed points. In this work, using redefinition of variables, we show that ELKO cosmology 
admits early-time stable fixed points. More interestingly, we show that ELKO cosmology 
admit two sets of attractor points corresponding to slow and fast-roll inflation. The fast-roll 
inflation attractor point is unqiue for ELKO as it is independent of the form of the potential. 
We also discuss the plausible choice of interaction terms in these two sets of attractor points 
and constraints on the coupling constant. 
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1 Introduction 

Inflation is currently one of the successful paradigm of the early universe[l, 2]. The success 
of inflation not only rests on solving the problems of the Friedman- Robertson- Walker model, 
and that it generates the primordial spectra of scalar (density) and tensor (gravitational 
waves) perturbations. The temperature variations of CMB as measured in WMAP and 
PLANCK — to a large extent — confirm that the these primordial density perturbations are 
generated due to rapid expansion of the quantum fluctuations in the early universe. 

Current CMB measurements can at-most provide two physical quantities during infla- 
tion: (i) inflaton (dynamical field that dominates during inflation) potential and (ii) the first 
derivative of the potential in the observable scales [3]. It is still unclear what are the prop- 
erties of the dynamical field that drive inflation. It is usually assumed that the field that 
dominates is a fundamental scalar field. While it is the simplest, recently, there has been a 
surge of activity to look at the plausibility that the inflaton has an internal structure [4-6] . 
One of the key results of the spinor condensate models compared to the standard inflationary 
models is the prediction of running of spectral index in the slow-roll limit that is consistent 
with the current CMB measurements [7]. 

Even if the observations provide us information about the nature of the dynamical 
field that drive inflation, the initial state of the field will never be known [8-10]. Einstein's 



- 1 - 



equations are non-linear, hence, it is important to know what range of initial conditions of 
inflation can plausibly lead to inflation. For instance, it is possible that if the initial velocities 
of the background field are large then this will stop inflation [11]. 

The knowledge of initial conditions can provide crucial information about the nature 
of the fields and their interactions with the known matter fields [11, 12]. For instance, it is 
usually assumed that the inflaton is a fundamental scalar field. However, we do not know the 
nature of the scalar fields or its interaction with other fields. Similarly, it is not clear what are 
the properties of ELKOs and how they interact with the other fields. If the observations do 
provide evidence that the inflation occured due to the one of these fields, the initial conditions 
of these fields will provide information about the nature of interactions with standard model 
particles. This in turn can be useful for model building which can be verified in high-energy 
experiments. 

The above issues are relevant and imperative to the current acceleration of the universe. 
It is unclear what dynamical fields drive the current accelerating universe. Even if the 
observations reveal the nature of the field, as in the case of inflation, it is not possible to 
know the initial condition that lead to current accelearation. This is referred to as cosmic 
coincidence problem. The constraints on the interaction of these fields with standard model 
particles will provide information about the initial condition that lead to acceleration. 

In this work we investigate the following questions: If dynamical field during inflation 
is a condensate of the non-standard spinor whether a large set of initial conditions lead to 
inflation. If it does, then, can it constraint the interactions between the spinor fields and 
the standard model particles. Recently a spin-half fermion called ELKO was proposed as a 
candidate of dark matter [13, 14]. Recently, in the literature several authors have shown 
that these can lead to accelerated expansion [7, 15-20]. 

In Ref. [21, 22], the authors could not find any stable fixed point with various kind of 
potentials. One of the draw-backs of their analyses is the choice of variables. Specifically, they 
have assumed <p an d V(cp) to be independent variables. In this work, we make a combination 
of <f>, H and V((p) and show that in these newly defined variables the dynamical equations 
have stable fixed points for a wide class of potentials and interactions between ELKO and 
matter. More interestingly, we show that ELKO cosmology admit two sets of attractor points 
corresponding to slow and fast-roll inflation. We also discuss the constraints on the coupling 
constants that will lead to early-time attractor behavior. 

The manuscript is organised as follows: In sections 2 and 3 we define the background 
equations and slow-roll parameters. In section 4 we define new sets of variables and rewrite 
the background equations in terms of these variables. A general analysis of fixed point 
behavior is presented in section 5 and section 6 contains the fixed point analysis for specific 
potentials. Finally, in section 7, we summarize our results and also discuss its implications. 

2 Background equations 

The new class of spinors named Elko were proposed as a candidate of dark matter in Ref. [13, 
14]. These spinors are eigen spinors of charge conjugation operator. These are non standard 
spinors(NSS) because unlike the standard spinors Elko spinors have mass dimension one and 
(CPTY = -I. As a consequence these spinors follow Klein-Gordon equation. Elko spinors 
are also called 'dark' spinor as its dominant interaction with standard model particles is via 
Higgs and gravity. Apart from being considered as inflaton these spinors have drawn lots of 
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interests in works of different authors [23-27]. The NSS(A) and its dual(A) can be defined as 

A = A = ml (2-1) 

where £ and £ are constant spinors with the property ££ = I, and <f){t) is a scalar which is 
time depepndent. The action can be written as [boeh] 



S = Id xy/—g 



\sT^{^u)X-V (AA 



(2-2) 



where V M is the covariant derivative and the round bracket in the subscript denotes sym- 
metrisation. Varying the action (2.2) the energy density and pressure can be written in terms 
as 

P^ = ^ 2 -^(^ 2 )'-^V-^), (2-3) 
where H is the Hubble parameter. The equation of motion can be written as: 

<£ + 3#0-jj# 2 + ^ = O, (2.4) 
3 Slow-roll parameters for ELKOs 

Due to the presence of the H 2 <fi 2 term in the density and pressure, one has be careful in 
defining the slow-roll parameters for the ELKO condensate. In this section, we give the 
expressions for the slow-roll parameters. From the expression of time-time component of 
Einstein's equation one can write the Friedman's equation for NSS as: 



3M 2 \2D 



V 



where D = \ — F=l — -f,* and V = X*. Taking the time derivative on the both the sides 

pi 

of equation (3.1) one can write the slow-roll parameter e as 

H 3 <j) 2 /D D , . 

e = ~TF? = ^^T, ~ + 777=: = £ can + a. (3.2) 

H 2 2(f)2/2D + V HD 

where a = -jjjj < 1 is an extra parameter which appears in the definition of the slow-roll 
parameter for ELKOs and e can is the slow-roll parameter defined for the canonical single 



scalar field, e can = § ■ ^ ~ |^ (when 4> 2 <C V). 



2 a jfl 

/2+V 

Substituting (4> 2 /2D + V) from (3.2) in (3.1) one can write 

2M 2 l H 2 D (e - a) . (3.3) 
Now in this case we define 5 = Taking the time derivative on both the sides of (3.3) we 



get, 



A = 5 = -e+*+ (£ 7 a) \ (3.4) 
H<f> 2 2H{e-a) V ; 
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The last term can be dropped as it is the time derivative of the slow-roll parameters. Therefore 
finally one can write the definition of 5 as: 

. a a 

5 = -e+ - = -e can - -. (3.5) 

A closer inspection of above expression immediately suggests that 6 is negative definite. For 
canonical scalar field, it is positive definite. This has an important effect for the spectra of 
scalar perturbations [7]. 



4 Dynamical equations 

The expressions of energy density and pressure (2.4) can be written in terms of newly defined 
variables X and V as following 

P4 = X + V, P4> = X-V, (4.1) 

where 



X = -tf - -{Htfy (4.2) 
V=^(H^)- + ^H 2 ^ + V(4>). (4.3) 

Physically, X refers to the kinetic energy of the condensate field and V refers to the potential 
energy of the condensate. Note that Hcjxf) acts like a friction term while Hcf) 2 acts as an anti- 
friction term. Elko dynamics crucial depends on which of these two terms dominate during 
inflation. 

Friedmann equation can be written as 

K 2 K 2 

H 2 = —Ptot = y (p<t> + Pra) , (4.4) 

where p m is the matter density and k 2 = 8ttG. Using equation (4.1) we can write the 
Friedmann equation(4.4) as 

x 2 + y 2 + v 2 = 1, (4.5) 



where x,y and v can be defined as x = V = and v = ■ Now, if we consider 

that the matter and dark energy are interacting only with themselves then the continuity 
equation 

p tot + 3H{ Ptot + p to t) = (4.6) 
can be written as two separate equations 

p <f) + 3H(p^+p^ = -Q, (4.7) 

p m + 3H{p m + Pm ) = Q, (4.8) 

where Q is the interaction term. In terms of the variables x, y, v equations (4.7,4.8) can be 
written respectively as 
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v '= ( e -^Av + ^. (4.10) 

Here ' is the derivative with respect to time divided by H, e = — ^ and A = ^, Q\ = 
To derive the above equations we have used the relation p m = (7 — 1) p m , where 7 takes the 
value either 1 or | depending on whether the universe is filled with cold matter or radiation 
respectively. Derivative of the variable y with respect to time give us 

+ (4 ' n > 

H can be written as 

K 2 

H = — \P4> + Pcf, + Pm + Pm] 

Therefore we have three dynamical equations (4.9), (4.10) and (4.11) with one constraint 
(4.5). It is important to contrast the above set of variables with those used earlier [21]. The 
two variables X and V are independent of each other. However, in Wei's analysis [21], the 
two variables y and u are not independent. 

In the rest of this work, we study the stability of fixed points with these equations. We 
show that ELKOs show a new set of fixed points that can not be identified directly with the 
canonical scalar field. 

5 Fixed points and stability analysis: General Analysis 

Fixed points are those points where the dynamical variables stop evolving, i.e., at fixed point 
(x, y, v) the time derivative of x, y and v are zero. At fixed points, dynamical equations (4.9, 
4.10, 4.11) can be written as: 

(e-3)x-Af _^ = 0j (5.1) 
H x x 

- ^7^ + ^ = 0, (5.2) 



e -+^).- = 0. (5.3) 
Eq. (5.3) leads to two set of fixed points: 

1. Case I: y = and e 7^ —jjj 

2. Case 11: y ^ and e = — 

In the rest of this section, we will consider the above two cases for general interaction term 
Q\. In the following section, we consider special cases for the interaction term and discuss 
the nature of fixed points. 
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5.1 Case I 



Substituting this value of y in equation (5.1) we get 

e = 3x 2 + 

x z 



(5.4) 



The above form of e gives crucial information about the class of interaction terms between the 
condensate and matter fields that can lead to attractor behavior. In particular, it immediately 
shows that Q\ oc x 2 may not lead to stable attractor points. Also it provides an upper bound 
on the coupling constant. We discuss these in the next section. 
General expression for e can be written as 



_H__3 

H 2 ~ 2 7 + 



2 1V 



(5.5) 



Therefore at fixed points e can be written as e = §7 + (3 — §7) x 2 . e is a positive quantity 
ensuring accelerated expansion of the universe. Finally one can write an important relation 
for e which will be used later 



2 7 



3) (l — x 2 ) . 



(5.6) 



Once we get the fixed points, we need to study the stability of the fixed point to ensure 
that the fixed points are actually giving us an attractor. If the fixed points are stable then 
we can have the attractor and finally we will be able to alleviate the 'cosmic coincidence' 
problem. To analyse the stability of these fixed points we perturb the system about the fixed 
point, x — > x + Sx and y — > y + Sy. Then we study the evolution of the perturbations. If 
we have a growing solution of the perturbations our fixed points are not stable, however if 
we find a decaying solution we can say that our fixed points are stable. Substituting these 
values of x and y in equation (4.9) and (4.11) we get the perturbed equations of x and y as 
follows: 



Sx' 



(e- 3) + (6 
ldQi 



3 7 )* 2 + ^-- 



Qi l dQi 

x dx 



Sx 



x dy 



Sy, 



(5.7) 



Sy' 



(5.« 



Here we have used Se = [(6 — 37) x] Sx and y = 0. Equations (5.7) and (5.8) can be written 
as 



Sy'J- {M) \Sy 



(5.9) 



where 



M 



(e --3) + (6-3 7 )^ + ^-I^ \ 
Two eigenvalues of the matrix M are 



Mi 



e + 



A 

2H 



(5.10) 
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M2 = (e ~ 3) + (6 - 3 7 ) x 2 + % - (5.11) 

X z X ox 

Stability around the fixed points depend upon the nature of the eigen values [i\ and /i2- 
When hi < 0, fi2 < the fixed points are stable and we can get an attractor solution. If 
fi\ > 0, fi2 > 0, the fixed points are unstable and we can not have any attractor. If one 
of them is positive and other one is negative, we get a saddle point which says that at one 
direction the fixed points are stable and at the other direction the fixed points are unstable. 
In Ref.[21] it was noted that equation-of-state parameter = y- > — 1 when cp 2 > 

I (H<ft 2 y. The dark energy will enter the phantom region (w^ < — 1) if 4> 2 < \ i^Hip 2 ^' . 
Therefore in the region > — 1 we always get X > 0. Now from Freidmann equation (4.4) 
we get 

H 2 = ^f(x + V + Pm ), (5.12) 

which implies that H 2 > ^pV". Finally taking logarithmic time derivative on both the sides 
of this inequality we get 

e + A <0 . (5.13) 

This means that [i\ is always negative as far as > — 1 is concerned. Therefore it is 
possible to have a stable fixed point if ^2 becomes negative for some interaction Q. In the 
next section, we analyse this for three types of interactions. 

5.2 Case II 

Substituting the value of e and constraint (4.5) in Eqs. (4.9) and (4.10), we get for (7 = 1): 

x ' = (e -_ 3)x _4y!_^i (5.14) 
H x x 

I V 



Substituting for A, we get, 



e = - A = 3x 2 + ^ 2 (5.16) 
Se = Qx5x + 3v5v (5-17) 



The perturbed equations about the fixed point are: 

„2 „,4 

Sx 



3-9x 2 - —v 2 - 3^ + 3^) Sx 12— - I5xv) 5v - 5(— ) (5.18) 

2 x z x z J \ x x J x 

Sv = 6 x v Sx + ^3x 2 + ^v 2 - Sv + S (Q^J (5.19) 

These attractor points are unique to ELKO cosmology regarding which we would like to 
stress the following points: 

1. The perturbed equations do not explicitly depend on the potential. Hence, these equa- 
tions can be realised for any potential provided e = —\/{2H) is satified. 

2. If x 1 and v — > 1 (or vice- versa), Eq. (5.16) implies that e > 1. It is easy to see that 
a(t) > implying that this corresponds to fast-roll inflation [28]. 
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6 Special cases of the interaction term 



In the previous section, we have obtained the condition for the existence of fixed point for 
general interaction. However, the analysis for a general interaction term is complicated. Here, 
for two cases, we take simple form of the interaction term and show explicitly the nature of 
the fixed points. 

6.1 Case I: Slow-roll 

6.1.1 Qi = (5v 2 x 

In this case the fixed point y is zero. The fixed point x and v can be found using the equation 
(5.1). Substituting y = in equation (5.1) and using v 2 = 1 — x 2 from equation (4.5) we can 
write 

(e-3)x-^ = 0. (6.1) 



Which gives us two solutions for x 



±1, % = - ? (6.2) 

( 3 - rn 



Now x = ±1 can not be a scaling solution because that will make our universe completely 
dark-energy dominated. Therefore the only possible solution is x = — i ^ 3 ^ which is negative 

as /3 and (3 — | 7 ) are both positive. 

Therefore in this case the eigenvalue (12 of the matrix M can be written as 



M2 



3 - | 7 J (1 - 3x 2 ) + 2/3s 



(6.3) 



Now substituting the solution of x from (6.2) in the above expression of \i2 one can get the 
following expression of \i2'- 

Prom the above expression of [12 it can be understood that when the first term dominates 
over the last term one can get ^2 < 0. Therefore the condition for having a stable fixed point 
for this kind of interaction is: 



Therefore in this kind of interaction the coupling constant can not be very large and e is 
always less than 1. 

6.1.2 Qi = /3v 2 x 2 

In this case equation (5.6) and (6.1) tell us that 

^7 - (1 - x 2 ) x - (3 (1 - x 2 ) x = 0. (6.6) 

So, the only solution of x = (0,±1). One can not take these solutions as fixed points as the 
universe will become purely matter dominated and dark-energy dominated respectively in 
those cases. Therefore we do not have any physical fixed points in this kind of interaction. 
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6.1.3 Qi = (3vx 2 

Following the similar method as described above using (6.1) for this kind of interaction one 
can find that at fixed point the only solution for x is: 



& 2 

?, (6-7) 

(3-§7f 

Here we have considered x / (0, ±1). Substituting the above expression of Q\ in the expres- 
sion of [i2 one can get 

x 2 

/u 2 = (6-3 7 )x 2 + /3— . (6.8) 



v 

Using the definition of v 2 = 1 — x 2 and the expression of x from (6.7) one can write the 
expression of [i2 in terms of the coupling [3 as: 

^ = V 3 "2 7 J-(3-|^)' (6 - 9) 

Therefore in this case the ^2 will be negative only when /3 > (3 — | 7 ) . However from (6.7) 
one can see that this condition will make x imaginary. Therefore we can not find a stable 
fixed point in this case. 

6.2 Case II: Fast roll 

6.2.1 Qi = (3v 2 x 

For this interaction, the perturbed equations of x and v are: 

Sx = [3 - 9x 2 - —v 2 - 3^tt + 3^}5x + [6- - 12— - Ibxv - 2f3v]5v. (6.10) 
2 x z x A xx 

1 9 3 

5v = [6xv + /3v]5x + [3x 2 + -v 2 - - + (3x]5v (6.11) 

The two eigen- values corresponding to the above set of equations are negative. Fig. (la) 
shows that for different initial conditions v — > 1 and x <C 1 is an attractor point. 

6.2.2 Qi = (3v 2 X 2 

For this interaction, the perturbed equations of x and v are: 



Sx = [3 - 9x 2 v 2 - 3^r + 3^r - (3v 2 ]5x + [6- - 12 15xv - 2Bvx]Sv. (6.12) 

2 x z x z xx 

Sv = [6xv + 2f3vx]5x + [3x 2 + \ y2 ~\ + /3x 2 ]5v (6.13) 

Here again, both the eigenvalues corresponding to the above set of equations are negative. 
The eigenvalues are negative for all ranges of /? for which x and v are real. Fig. (lb) shows 
that for different initial conditions v — > 1 and x ^ 1 is an attractor point. 
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6.2.3 Qi = f3vx 2 

For this interaction, the perturbed equations of x and v are: 

Sx = [3 - 9x 2 - —v 2 - 3^r + 3^r - M8x + [6- - 12— - l$xv - /3x]5v. (6.14) 
2 x z x z xx 

9 3 

tfv = [Q XV + 2/3x]fo + [3x 2 + -u 2 - -]<5u (6.15) 

Here again, both the eigenvalues corresponding to the above set of equations are negative 
for all valus of f3 where x and v are real. Fig. (lc) shows that for different initial conditions 
v — >■ 1 and x <C 1 is an attractor point. 

7 Conclusion 

As it is known, it is not possible to know precisely the initial condition of the field that drives 
inflation. Hence, any model of inflation need to show that attractor points exists in the space 
of the matter field variables. In all previous work, it was not possible to show explicitly that 
the Elko condensate does indeed lead to late-time attractor behavior. 

In this work we have shown that rewritting the background field equations interms of 
new variables X and V , one could show the existence of late-time attractor. Interestingly, 
we notice that the two set of attractor points exist for ELKO condensates. In case I, which 
is similar to canonical scalar field inflation [11], the attractors are realised only when e < 1. 
In case II, which is unique for ELKO cosmology, the attractors are realised only when e > 1 
and they exist independent of the form of the potential. It will be interesting to repeat 
Starobinsky's analysis [28] for this case and see whether the spectrum of perturbations are 
nearly scale- invariant. 

In the phantom region (wj, < —1) obtaining stable fixed point may not be possible as 
in that case X < may not allow H 2 > V, in other words in phantom region fj,i < 
may not be possible always. In this work the variables we have worked with are not a 
simple transformation of the variables chosen by Ref.[21]. The stability may be achieved by 
redefining the potential and the kinetic part. This is currently work under progress. 
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